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The improved quantum scheduling algorithm proposed by Grover 
has been generalized using the generalized quantum search algorithm, 
in which a unitary operator replaces the Walsh-Hadamard transform, 
and 7r/2 phase rotations replace the selective inversions, in order to 
make the quantum scheduling algorithm suitable for more cases. Our 
scheme is realized on a nuclear magnetic resonance (NMR) quantum 
computer. Experimental results show a good agreement between theory 
and experiment. 
PACS number(s):03.67 

1 Introduction 

The scheduling algorithm solves the intersection problem, i.e., the problem 
is to find common elements in two sets. Let Alice and Bob be two distant 
parties who wish to collaborate on a common task. They each have a schedule 
listing N slots of time. Their schedules can be represented as two strings if 
Alice and Bob each have N classical or quantum bits. For each person, a bit 
in state 1 denotes a slot that is available, and a qubit in state denotes a slot 
that is not available. The intersection problem becomes to find the common 
Is in the two strings. Alice and Bob need to exchange information to find 
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the common Is. The problem is how to reduce the exchanging information 

m 121 ia 

If Alice and Bob exchange classical bits, they will need to exchange O(N) 
bits. H. Buhrman et al found that Alice and Bob could find the common Is 
by exchanging 0(VN \og 2 N) qubits, if they exploited the quantum paral- 
lelism jHIHj- Grover proposed an improved algorithm in theory for the special 
case in which the string of either Alice or Bob has few Is 0. Grover's the- 
oretical scheme can be described as follows. Assume that there is a log 2 N 
qubit register, by setting N = 2 n , where n is an integer. The register can be 
transmitted between Alice and Bob. Assume that Alice has eiV Is in her N 
bit string, where eN <C N. For convenience, we assume that Bob also has 
eN Is in his N bit string, and there is a single common 1 in the two strings. 
Alice encodes the N slots in the register, by applying the Walsh- Hadamard 
transform (denoted as W) to |0 >, where |0 > denotes that all qubits in 
the register lie in state |0 >. After she repeats WIqWIa tu times, where 
m = irJlJe/A [7j, the register lies in the superposition corresponding to Al- 
ice's available slots. I a denotes the selective inversion (or suitable inversion) 
for the basis states corresponding to Alice's available slots. Iq denotes the 
selective inversion for |0 >. The composite transformation used by Alice is 
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denoted by G, represented by G = (W IqW I A) m W ' . According to the gener- 
alized quantum search algorithm jH], (GIqG^Ib^G transforms |0 > to the 
state corresponding to the available slot of the both of Alice and Bob, where 
G^ 1 denotes the inversion of G, and Ib denotes the selective inversion for 
the basis states corresponding to Bob's available slots. 

Alice can carry out W, Ia, and Jq. When she needs Ib, she must send the 
register to Bob, who applies Ib to the register, and then returns it to Alice. 
One should notes that the overall state of the register is unaltered during the 
course of sending 0. Obviously, the number of times the register needs to 
be sent to Bob is equal to the number of Ib operations. 

In the improved scheduling algorithm proposed by Grover, he used the 
original quantum search algorithm, which is not efficient or invalid at all in 
some cases The limits on the search algorithm impose restrictions on 
the quantum scheduling algorithm. In this paper, we generalize Grover's 
improved algorithm by replacing W by a proper unitary operator U, and 
replacing the selective inversions by tt/2 phase rotations, in order to make 
the quantum scheduling algorithm suitable for more cases. Correspondingly, 
Jq is also replaced by a n/2 phase rotation for |0 >. We implement our 
scheme using nuclear magnetic resonance (NMR). 
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2 Generalizing the quantum scheduling algo- 



Our experiments use a sample of carbon-13 labelled chloroform dissolved in 
d6-acetone. Data are taken at room temperature with a Bruker DRX 500 
MHz spectrometer. The resonance frequencies V\ = 125.76 MHz for 13 C, and 
U2 = 500.13 MHz for 1 H. The coupling constant J is measured to be 215 
Hz. The Hamiltonian of this system is represented as ^U] 



by setting H = 1, where J* (k = 1, 2) are the matrices for ^-component of the 
angular momentum of the spins. The evolution caused by a radio-frequency 
(rf) pulse on resonance along x or — y axis is represented as R^.(^p) = e tipI ^ or 



is denoted by or [— ip\y. The coupled-spin evolution is denoted as 
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Ry(—(f) = e t(pI y , with k specifying the affected spin. The pulse used above 
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where r is evolution time. The pseudo-pure state 
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100 > = 
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is prepared by using spatial averaging 1 1 j , where |0 > denotes up spin state. 
The basis states are arrayed as 1 00 >, |01 >, |10 >, |11 >. In experiments, 
we exploit rf and gradient pulse sequence to transform the system from the 
equilibrium state to the state po = + If/2 + l\l 2 z -, where po denotes the 
deviation density matrix equivalent to 1 00 > [P2] . 

The Walsh-Hadamard transform is replaced by U = Ry(ir/2)Ry(—7i/2), 
which is represented as 
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Alice and Bob each have a four bit string, in which the positions of slots 
are encoded by a two qubit register. The basis states 1 00 >, |01 >, |10 >, 
1 11 > correspond to the first, second, third, and fourth slots, respectively. 
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Assume that Alice's available slots are the first and the second ones, which 
correspond to 1 00 >, and |01 >, respectively. I a is replaced by a n/2 phase 
rotation jT^j, represented as 
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Considering the phase matching condition [T3] , Iq is replaced by 
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We assume that there are three cases for Bob's available slots: 1) The 
first and the fourth slots are available; 2) The second and the third slots are 
available; and 3) The first and the second slots are available. For the three 
cases, Ib is replaced by 
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respectively. By denning the composite operators G i2 = — C/ _1 /| 12 C/, 

Qi = — G\2^ G12 IbiaGi2, Q2 = — G\2^ C 1 2 1 /J 2 3 < ^i2) and Q 12 = — G^I^ G^I^G^, 

we obtain 



G 12 |00 >= e- J7r/4 (|00 > +|01 >)/v / 2, 



(10) 
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QijOO >= — iJOO >, 



(11) 



g 2 |oo >= — i|oi >, 



;i2) 



g 12 |00 >= e"^ /4 (|00 > +|01 >)/v / 2 



(13) 



Eq. fjllj) . for example, shows that the register is transmitted between Alice 
and Bob only one time, and the common slot, which corresponds to 1 00 >, 
is obtained. Similarly, if Alice's available slots are the third and fourth ones, 
we replace J| 12 by 
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and obtain 



G34IOO >= - e - i7r/4 (|10 > +|11 >)/v / 2 



(15) 
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where G 34 = —UI^U~ l IX 3A U. We also assume that there are three cases for 
Bob: 1) The first and the fourth slots are available; 2) The second and the 
third slots are available; and 3) The third and the fourth slots are available. 
By defining Q 4 = -G M I-2 G 3 il£ u G u , Qs = -G M I^ G 3 il£ 23 G u , and Q u = 
-G 3 J^G 3 il^ u G 3 4, we obtain 



where 



g 4 |oo >= >, 



(16) 



g 3 |oo >= zjio >, 



(17) 



g 34 |00 >= -e- m/4 (|10 > +|11 >)/v / 2, 



(18) 
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The overall phases before wave functions can be ignored. 
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3 Experimental procedure and results 

The coupled-spin evolution described as Eq.© is realized by pulse sequence 
t/2 — [ti]1 ,2 — t/2— [— 7r]^ 2 [IB], where [ti] x 2 denotes a nonselective pulse (hard 
pulse), and the symbol t/2 denotes the evolution caused by the magnetic 
field for t/2 when pulses are switched off. [n] x 2 pulses are applied in pairs 
each of which take opposite phases in order to reduce the error accumulation 
causes by imperfect calibration of 7r-pulses [T7]. U is realized by [tt/2]* — 
[-tt/2] 2 . /f 12 = i2j(7r/2)i£(7r/2).Rj(-7r/2) (up to an irrelevant overall phase 
factor), and is realized by [—n/2] y — [tt/2]]. — [n/2] y , noting that the pulses 
are applied from left to right [TH]. IJ 14 = [7/2 J], JJ 23 = [1/2J], and 7J 12 = 
I\xi- By modifying the pulse sequences in Ref. (THl Ej, we find = 
Rl> 2 (n / 2) Rl> 2 (n / 4) R y > 2 (-it /2)[15/ 4 J}. By optimizing the pulse sequences, 
we obtain 

G 12 = -Rl(7i/2)Rl(n/2)R 1 x 2 (n/A)R 1 /(-7v/2)[15/AJ}Rl(7v/2), and 

Gn = -i?K-V2)[l/4J]^ 2 (7r/2)^ 2 (-7r/4)^(-rr/2)i?i(-7r/2). 

Similarly, ij 34 = R 1 y (ir/2)R 1 x (-7r/2)R 1 y (-n/2), and lj 34 = /f 34 . G u can 
be optimized as G 34 = -R 1 y (7r/2)R 1 y (iT/2)R 1 x 2 (7T/A)R^ 2 (-TT/2)[15/AJ]R 1 x (-7i/2) 

The experimental results are represented by the deviation density matri- 
ces reconstructed from the spectra recorded through the readout pulses by 
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the technique of state tomography (T3]. We first prepare pseudo- pure state 
1 00 > . The generalized quantum scheduling algorithm starts with 1 00 >. Qi, 
Q 2 , and Q X2 transform |00 > into |00 >, |01 >, and (|00 > +|01 >)/v / 2, re- 
spectively. Figs.l(a-c) show the experimentally measured deviation density 
matrices after Qi, Q 2 , and Qi 2 are applied to 1 00 >, respectively Figs.l(d- 
f) show the theoretically expected matrices 1 00 >< 00 1, 1 1 >< 01 1, and 
(| 00 > +|01 >)(< 00 1+ < 01|)/2, respectively In contrast with the theo- 
retical expectation, for the values of the theoretical non-zero elements, the 
relative errors of the experimental values in Figs.l (a) and (b) are less than 
12%, and the relative errors in Fig.l (c) are less than 22%. The relative 
errors increase in Fig. 1(c) because the theoretical values are only half of 
those in Figs.l (a) and (b). The other small elements in Figs.l(a-c) are less 
than 30%. The errors mainly result from the imperfection of pulses, effect of 
decoherence and inhomogeneity of magnetic field. 

4 Conclusion 

We have realized the generalized quantum scheduling algorithm on a two 
qubit NMR quantum computer, and obtained nontrivial results. For each 



12 



case discussed above, Alice and Bob each have two available slots in their 
schedules either of which consists of four slots. The original scheduling al- 
gorithm does not work in this case, because the original search algorithm 
is invalid in the case that the number of the marked states is equal to the 
number of unmarked states. However, using the generalized scheduling algo- 
rithm, Alice and Bob can find the common slots by exchanging the register 
only one time. 
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Figure Captions 

1. Experimentally measured deviation density matrices shown as Figs. 1(a)- 
(c) and theoretically expected density matrices shown as Figs.l(d)-(f) 
after the completion of the quantum scheduling algorithm, correspond- 
ing to operations Qi, Q2, and Q12, respectively. Only the real compo- 
nent is plotted. The imaginary portion, which is theoretically zero, is 
found to contribute less than 15% to the experimental results. 

[Figure 1 about here.] 
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